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We review one dimensional matrix theory and its variations, collective field theory and quantum 
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I. INTRODUCTION 


The idea of large-A^ limit was implicitly proposed in early 1950s, for different models 
in condensed matter physics. See [l| for the early history and development. In 1974 ’t 
Hooft applied it to QCD j^, and reorganize the expansion of the theory in terms of 1/iV 
and A = Qym^- The leading diagrams in 1/N can be drawn only on a plane/sphere, 
and termed planar diagrams. The large- iV limit with such a reorganized expansion is thus 
known as planar limit. The 1/A^-suppressed terms can only be drawn on surfaces with higher 
genus, which indicates some relation with the topological expansion of the string scattering 
amplitudes. For N = 4 Super Yang-Mills theory, the planar diagrams are described by the 
non-interacting IIB string theory on 5 dimensional anti-de Sitter spacetime multiplied by a 
5 sphere j^. They are given as an expansion at large’t Hooft coupling A, with the leading 
contribution generated explicitly by the corresponding supergravity. Following ’t Hooft, 
the planar expansion is extended to the scalar field theories, where the scalar is uplifted 
to a matrix field jd]. In one dimension, the planar limit of such a matrix theory is solved 
by a non-interacting fermionic system, with the fermion positions representing the matrix 
eigenvalues. To describe realistic many-body system at large- iV, one could generalize the 


matrix theory to a theory of the density function of the particles, the collective field 


a. 


ensity operator, which 
8|, it is shown that in 


Alternatively, one could use the phase space representation of the d_( 
is called W^o coadjoint orbit or coherent state approach [g, [TJ. In 
general the coadjoint orbit/coherent state approach captures the classical dynamics in the 
large-Y limit. 

Despite the early development of large-Y limit in condensed matter physics, the planar 
expansion is not widely employed. In 2009, Ma and Yang analyzed the Y-dependence 
of the groundstate energy of ID two-component fermionic system, with both a ^-function 
interaction and a harmonic conhning potential . They found that the energy scales as Y^ 
when the interaction is repulsive, and —g‘^N when the interaction strength g —)■ —oo. They 
conjectured that the re-scaled energy E/N'^ should be a smooth function of g/y/N- One 
may immediately recognize this as a planar limit, of a properly generalized matrix theory. 
In this paper we will show how this could be indeed realized. 

We organize the paper in the following way. In the next section we review the ID matrix 
theory, its planar limit, and the solution through non-interacting fermions. In section HI 
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the collective field theory is briefiy shown, and applied to ID interacting bosonic system. 
Then in section IV we nse the phase space description to describe ID interacting fermionic 
system. We take the 5-function interaction as an example, and compare the results with 
those from other approaches when available. In the last section a short summary is given. 


II. MATRIX THEORY, PLANAR DIAGRAMS AND ID FREE FERMIONS 


First we review the relation between the ID hermite matrix theory and the free fermion 
system jd]. We will repeat the detailed procedure, in order to show the relation to the 
fermionic system clearly. For the 93 ^ theory, the matrix Lagrangian is given as 


L = ti{dtMdtM^) + tr(MMt) + (1) 

Following the proposal of’t Hooft , the propagators of the matrix field are represented with 
double lines. Consider a general connected vacuum diagram, which contains P propagators 
, V vertices and I closed loops of internal index. Viewing the internal loop as a surface, one 
obtains the Euler relation 

V - P + I = 2-2H. (2) 

Here H is the number of holes of the Feynman diagram. Taking into account 2P = 4V, one 
can re-arrange the factors to obtain the weight coefficient for the diagram 




( 3 ) 


Therefore if we take oc 1/N, only planar diagrams with H = 0 remain in the large-iV 
limit. In other words, one only needs to keep those diagrams that can be drawn on a plane 
or sphere. Explicitly, one may write 

= (4) 

where Eolg^) is the sum of all connected vacuum diagram on a planar surface. Notice that 
the conclusion will not be changed even if the interaction is of power p 7 ^ 4. In such a case 
the corresponding coupling constant scales as gp oc In particular, for the quadratic 

term the coefficient does not change with N. It simply reflects the difference between the 
propagator and the interaction. 


lim / 

N—yoo J 
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In general, it is still not easy to snm all the planar diagrams. In the present ID case, it 
turns out to be equivalent to hnding the ground state energy of the corresponding Hamilto- 
mian 

^ ^ ^ 2 ^ VaReM2 + dluiMl) 

V = ^trM^ + 5(4trM^. (5) 

That is, we only need to solve the equation 

Hil^ = N‘^E{g,)'il^, (6) 


and hnd the ground state wave function. The wave function should be symmetric under the 
U(N) rotation M —)■ UMU~^, as required in the matrix theory. Alternatively, the ground 
state energy can be obtained from minimization of the action within the invariant wave 
function conhguration 


Eig^) 


lim —- min 

N^oo p 




( 7 ) 


Therefore it indeed corresponds to the sum of the connected vacuum diagram (j4]). The 
integration over the angle part U can be done trivially, leaving the integration over the 
eigenvalues Xi of M. So the ground state energy becomes 


Ui ^XiUi<jixi - Xj) 


E{g 4 ) = lim ^ min ■ 

N^oo iy ^ Ip 




n* dxi n* 


<j 


.Xi 




Dehning 


i){xi, ...,Xn) 



one obtains the Schrodinger equation 


il^Xi, ...,Xn), 




( 8 ) 

( 9 ) 


( 10 ) 


From the dehnition ([9]) we see that the above equation describes the emotion of N fermions 
in ID. They are conhned in a central potential, but do not interact with each other. In 
the large-A^ limit, they behave semi-classically and simply £11 all the states below the Fermi 
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energy level ep- Integrating out the momentum, the energy and the fermion number are 
then expressed as 


N^E{g,) 

N 



— — 25(40;^]^/^ 9 { 2 eF — — 

6{2eF — x^ — 2g4X^) 


2giX^) 


( 11 ) 


Clearly at lage N, ep ^ N, and the range of x, as well as the density distribution of x, scales 
as y/N. Re-scaling the parameters as ep = Ne and x = y/Nu, the above relations can be 
simplihed 


E{\) 

1 


e — [ — [2e — 6{2e — — 2Xu^) 

J Svr 

J — [2e — 6{2e — — 2 Am^), 


( 12 ) 


where A = g 4 :N. The density in the new variable u now reads 

p{u) = —\^e — V? — 2 \u^Y^‘^ ^(2^ — V? — 2\u^). (13) 

TT 

Taking (74 = 0 one recovers the famous semi-circle law of Wigner , after redehning u to 
absorb £. When g^ is non-vanishing the semi-circle is deformed. One may set the quadratic 
term to zero instead, then the distribution function is completely determined by the quartic 
term. Nevertheless, scaling of various quantities with N is unchanged. 

How to generalize the above discussion to bosons? How to introduce the interactions 
between the particles? It is not so obvious in the matrix formalism, although some progress 
could still be made. In the next two sections we will show how to deal with them in the 
so-called collective held theory and the quantum phase space description. Before plunging 
into details, we can hrst sketch the iV-scaling when the interactions are present. Since the 
coefficient of quadratic term (mass) does not change in the planar limit, the range of x and 
the density function 4>{x) in x are always of order yfN. Integration (f){x) over x then gives 
the correct particle number N. As a result, for a system with the ^-function interaction 


AH = gJ2^i^i- 


X 


3)1 


(14) 


i<j 


the planar limit can be achieved when g/y/N is kept hxed. This is exactly the hypothesis 
made in {q] and further conhrmed in H, 12|. With the Wscaling of the coupling con¬ 
stant taking into account, we could also recover the large-behavior obtained for a power 


potential trap 


13|- 
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III. COLLECTIVE FIELD THEORY AND ID INTERACTING BOSONS 


First we ignore the problem of the statistics, and try to reexpress the 5-fnnction interac¬ 
tion in a matrix form. The answer is almost immediate 

tr(e-'“') tr(e'“'), (15) 

i<j ^ 

However, the factor from the matrix integration measnre ([8]) will force the 

contact two-body interaction to vanish. Therefore, it is convenient to leave the matrix 

formalism and deal directly with the collective held |5| 

4>{x) = / = ^6{x - Xi). (16) 

It describes the density of the particles at a hxed position x. Now we want to nse (j){x) to 
describe the system, instead of the positions xi,...,xn- Explicitly, one could restrict x in 
some hnite interval —Lj^ < x < L/2, and start with the Fourier modes (pk of 4>{x) 

(j)k = ^ j dxe~^’''^(p{x) = itr(e"‘''^), (17) 

with k = 21:711 L and n takes integer numbers. Clearly 0(a:) contains more than enough 
degrees of freedom than the positions, even for hnite L. That means the Fouries modes 
(pk may not be independent. However, at least in the large-iV limit this will not cause any 
problem. 

Consider a bosonic system with a general Hamiltonian 

N N N 

7/=-^p2^-^n(xi,Xj) + ^fr(a;i). (18) 

2=1 27 ^J 2 = 1 

We need to make the coordinate transformation from Xi to (pk- In other words, we consider 
now the many-body wave function pilxi] as a composite function '0[0fc(xj)], and rewrite H in 
terms of (pk- In making the variable change one has to take care with the Jacobian. Taking 
L —)■ cx) in the end, the expression in terms of (pk can be compactly expressed through the 
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B 


field (j){x). After a sightly long derivation one obtains the following Hamiltonian 

^2 


= J dx 

Ay = i 

Vo = / 


^d^TT(j)d^TT + + AH + Ho 

^^dx-Uo.m 


(19) 


Hy) 

(x-y) 


dy 


dx 


(/ii? — H(a;) + v{x, a;))0(a;) 


dx 


1 

2 


di / da ij>(x)v(x, y)ij>(y). 


( 20 ) 


ensure 


where 7r(x) = is the conjugate held of 0(x), /Xi? is a Lagrangian multiplier to 

the constraint 

/ dx 0(x) = N. (21) 


Notice that the hrst term in AH is the well-known Weizsacker term 


15| . introduced to 


better describe the kinetic term in the Thomas-Fermi model. With the previous A-scaling 
arguments, one can easily hnd that AH is suppressed in the planar limit. This will be 
manifest if we rescale the density function and the coordinate as 


X = \/]Vm, 0(x) = '/Np{u), 


( 22 ) 


and redehne the coefficients in the potential H(x) and the interaction v{x,y) accordingly. 
The re-scaled variables u and p will be independent of N. Only an overall factor remains 
in the Hamiltonian, as expected from the discussion in the previous section. Since the held 
behaves classically in the large-A limit, we can determine the density function by minimizing 
the potential. Notice that although the kinetic term in the Hamiltonian is of leading order 
A^, the contribution can be neglected when we consider the density to be stable and not 
evolving with time. 

Let us apply the above formalism to an explicit example. We take the potential to be 
harmonic, and the interaction of the ^-function form 


(23) 


V(x) = -tx^x^, v(x,y) = gS(x - y). 

When the harmonic potential is absent, this gives the Lieb-Liniger model 
solved analytically through the Bethe ansatz . In the planar limit u: remains unchanged, 
and g behaves as \/A. Therefore the relevant coupling will be a = g/y/N. So the conjecture 


16| , which could be 


made in 


9| is just the consequence of the planar limit. An equivalent parameter has been 
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previously used in ly, [iSj]. Actually one may identify the large-limit in 1^ as the planar 
limit, with the density dehned there scaling as '/N. An immediate consequence of the 
above scaling fl2^ is, in each small interval dx with nonzero density, there will be as many 


as 0{\/N) particles. One could apply a localized version of the 
hnd the local energy density and Gibbs energy density llMllll. 


redholm equation 


16|, and 


18| . These results can then 


be used as input to solve the remaining hydrodynamic equations in a potential. 
With the new variable and held, the leading order potential reads 


Veff = 



^,.22 

/ 1Y'’ (“) - 



a 


pM + !-d“- 


(24) 


where 


jj,p = flp -\- Ci6{u — ti). 


(26) 


The symbol 5{u — u) could be thought of as the diagonal element of the identity operator 
in the coordinate representation. 


6{u — u) = {u\f\u). 


(26) 


As argued in 19|, to apply the Coleman-Luther-Mandelstam bosonization to the non- 


relativistic fermion, one has to extend the fermion sea to inhnity. The divergent a6{0)/N 
term for any non-vanishing a, though sub-leading in N, may provide a natural explanation 
for such an extension. Therefore as long as a is non-vanishing, the system will exhibit some 
kind of fermionic behavior, and the Thomas-Fermi approximation could be used 20j. As 
a —)■ 0, one would expect a phase transition to a true bosonic phase. 

Doing the functional derivative with respect to p{u), one obtains the corresponding equa¬ 
tion 


TT^ 1 

' 2 / \ / \ "^22 
— p [U) + ap[U) = Sa — -00 u , 


(27) 


This is essentially an explicit hydrodynamic equation [1^, |2^ . Despite the simple form of 
the above equation, it is not quite easy to solve due to the implicit dependence of £ and the 
range of u on a. This is hidden in the normalization 


J du p{u) = 1, (28) 

which is inferred from fl?ID . Moreover, both the limits a —)■ 0 and a —)■ cx) are not smooth 
and correspond to a phase transition. When a —)■ cx), the two a dependent terms are singled 
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out 


K// = 


a 


a 


p (w) - p(“) 


Functional differential to p{u) leads to 


P(w) = 

In other words, the corresponding density operator p satishes 


P = Jp, (NpY = Np. 


(29) 


(30) 


(31) 


In the next section we will show that this condition specify the fermionic nature of the 
system. When such a constraint is satished, the interaction term vanishes identically. This 
means that an inhnitely strong 5-function interaction is equivalent to a fermionic Pauli 
exclusive force, as hrst shown in 2l|]. The system in such a limit is called Tonks-Girardeau 


gas 


2l|, l22|. fl3T|) does not £x the function form of p{u), which is determined by the distance 
of successive particle positions. One then has to go to the sub-leading term in a. Repeating 
the procedure, one hnds 


p{u) = -\/2eoo - 

TT 


(32) 


Thus one recovers the semi-circle law for free fermions trapped in the harmonic potential (fl^ . 
The parameter Eoo will be hxed by the normalization fl25]l . giving The groundstate 

energy is easily obtained from by performing the integration in fl2T|) . with the result 


Eoo — N Eoo + Vgff — 


(33) 


This is simply the ground state energy of N spinless fermions in a harmonic potential, as it 


should be 


23j. One may easily generalize the above discussion to the original Lieb-Liniger 


model, with no conhning potential. At inhnite interaction, one hnds £00 = V) giving 


^00 = ^ N \ 
0 


(34) 


This is indeed the ground state energy for N spinless fermions at large N 211. 

The free fermion nature of the collective held theory at inhnite contact interaction allows 


for a much simpler derivation of the Hamiltonian (IT^ . 


one-particle density matrix element of the fermion held 


24 


'he collective held is then the 


26| 


(f){x) = T^(a;)\l'(a;). 


(35) 
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Through bosonization for the fermion field T, one could easily obtain the leading two terms 
in flT^ [w, 27|. As argued before, the divergent a6{0)/N term may be necessary to apply 
such a bosonization for non-relativistic fermion. The bosonization also suggests a 2D string 


theory representation of the problem 


27| |. Alternatively, one could deduce the dynamics 


of T directly in the large-A^ limit. Collective evolution of \h(x) is simply reflected in the 


equation 


mo 

ion 

24 

26 


26| . Translating the energy of the fermion liquid back with the above formula. 


one immediately obtains the leading two terms in flT^ . Since these two terms come from 
the kinetic part of flTS]) . they remains of the same form away from the free fermion point. 
Likewise, giving fll9p together with the constraint fl3ip . one could retain the Euler’s equation 


for an ideal fluid of free fermions with the following identification 


19| 


V = dxTT, 


P = \4>\x). 


(36) 


Here v is the fluid velocity, and P is the pressure. From this comparison one recognizes 
the close relation to the hydrodynamic approach, which is usually derived from the Gross- 
Pitaevskii equation 20|. Our previous approximation by neglecting the kinetic term corre¬ 
sponds to the static case there, with u = 0. And the conjugate field n^x) should be identified 
with the phase operator 6. 

When a —)■ 0, it has long been known that the density function is singular Q . Later we 
will show that the density goes to a 5-function. If a is small enough, the first term in the 
equation can be treated as a large constant with little variation. After absorbing it into the 
chemical potential e, one obtains the parabolic form density 


p{u) = - 
a 


Eq - 

“ 2 


(37) 


This is exactly the result obtained with the Thomas-Fermi approximation 2^ . See also the 
hydrodynamic derivation in and the numerical confirmation in 23]. The parameter Eq 
can be determined from the normalization fl28l) to be 


^0 — 


3aa;\ 

471 ; 


(38) 


The corresponding energy of the ground state is 

, 3 , 3 / 3uj \ 


2/3 


a 


2/3 ^2 


(39) 
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This is exactly the leading term found in 23 


together with a localized Fredholm equation 


obtained in the hydrodynamic approach 


11 


18| . The present framework therefore 


provides a natural basis for the validation of such a localized approach. When the harmonic 
potential is turned off, the chemical potential turns out to be £o = ctj resulting 




(40) 


Again without using the Fredholm equation, which results from the Bethe ansatz [17|] at 
large iV, we recover the leading term in the ground-state energy Q However, similar 
as the Fredholm equation, it seems not easy to go beyond the leading order analytically, 
especially at small a. 

Now let us take the a —)■ 0 limit. It is easy to find from flTTI) p(0) ~ while 

the maximum value for nonzero density, Uq, decreases as Taking into account of the 

normalization fl28p . one concludes that as a ^ 0, 


p{u) 5{u). 


This corresponds to a transition to the Bose-Einstein condensation (BEC) 
words, the BEC phase shrinks to a single point a = 0 in the planar limit 


(41) 


^ase. In other 


20 |, 


23 


28j. We 


will get the same conclusion if we start from N harmonic oscillators with no interaction, and 
then take the large-limit. In this case the effective potential is dominated by the limear 
terms 

1 r ^ ^ r . 

(42) 


I/O = - 

e// c 


(dx^) 1 2 f 

(P ^ J 


Due to the absence of the interaction, the chemical potential pp vanishes and Bose-Einstein 
condensation occurs. The ground state density function can be directly obtained from 
minimizing and is given by the Gaussian function 


M^) = N 


TT 


(43) 


Such a derivation is similar to the linearized Gross-Pitaevshii equation 20|]. It is not difficult 
to check that when iV —)• cx), 

(j)o{x)/'/N ^ 6 {x/'/N), (44) 


^ The normalization for the energy here differs by a factor of two from 161. 
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just as in eq. gill. As a result, the N'^ part of the ground-state energy fH2|) vanishes and 
the non-trivial contribution appears at 0{N). 

For intermediate a, numerical technique is needed to hnd the exact density function. 
When the density p{u) is obtained, we can again integrate flT^ to get the energy. Notice 


that we never deal with the momentum distribution, as studied in ll|, [ly, l23|. It will be 


interesting to start from a momentum density function p{k), and express the Hamiltonian 
completely in terms of p{k). One would expect that minimizing the Hamiltonian gives rise 
to some generalized Fredholm equation for p{k) [h, 2, However, it seems not easy to 
express the (5-function interaction in momentum picture directly. Perhaps the introduction 
of the phase space density is needed, as shown in the next section. 


IV. PHASE SPACE DESCRIPTION AND ID INTERACTING FERMIONS 

Now we try to extend the above formalism to the fermion system. As shown in 
in the inhnite interaction limit the density function (^(x) be expressed through the fermion 
held \b(x). For interaction fermions it is natural to consider the following operator 




(45) 


where |4/(t)) is the single-particle state vector. In the coordinate basis, one hnds 

^{x,y) = ^{x,t)¥{y,t). (46) 

Using the anti-commutation relation, it can be shown that <F(f) satishes the following relation 

<|(t)2 = (1 +v)<l>(t). (47) 

Employing the second quantization relation for fermions, the proper one-particle density 
operator could be dehned as jo] 

(^ = i - I), (48) 

which in coordinate representation reads 


(t){x,y) = 4'^(x,f)4'(|/,t). 


From the above expression it acquires the correct normalization 


(49) 


ti 6 = N. 


(50) 
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With fH7|) one can check that it satishes also 


(51) 


which represents the fermionic structure of the system. In the previous section, the operator 
Np in the strong-interacting limit satishes exactly the same constraint fl31l) . The free fermion 
system can be conveniently described through the coadjoint orbit/coherent state of the 
IToo algebra j^, which is the algebra of differential operators in the single-particle Hilbert 


space 


Roughly speaking, under the action of the Woo transformation, the operators 


follow a coadjoint orbit, and the states become co^ 
relation with the two-dimensional string theory 


lerent ones. Such a description has a close 


24 


26l |. For non-interacting fermions in a 


central potential, the action can be compactly written as 


S'[0] =i / dsdttT{^[dt(p,ds 


df tr(0h). 


(62) 


where h is the single-particle Hamilton operator, with the element 

(x\h{t)\y) = h{x,y,t) = hsj - V{x))S{x-y). 


(53) 


Varying cp along the adjoint orbit with (l5T|) and fl50|) preserved, one gets the equation of 
motion 

idtp + [h, 0] = 0. (54) 

In the time-independent limit, the equation can be easily solved. Using the phase space 
representation, the equation becomes 


[h(p, q),<t>{p,q)}MB = 0, 


(55) 


where MB denotes the Moyal bracket 30|]. Therefore (p{p,q) should depend on the phase 


variables only through h{p,q). Moreover, in the large-V limit the constraint (IHT]) simply 
requires 

(j){p,q)^ = (j){p,q). (56) 

Taking all these into account the ground state density is given by 


q) = ff((F-h(p,q)), 


(67) 


with ejp the Fermi energy determined by the conservation of particle number. 
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It will be interesting to extend the discussion to the interacting case. In particular we 
want to discussion the case with a h-function interaction as in the bosonic case. Due to 
the Pauli exclusion principle, such a interaction makes no effects among spinless fermions. 


Therefore one has to introduce the spin degrees of freec 
system with a ^-function interaction is first studied in 


om. The two-component fermion 


3l| and [3^, and thus called the 


Gaudin-Yang model. Recently the exact solution when the interaction strength g ^ oo is 


obtained in 


33l | . A nice review of the historic development of the model and the experimental 


progress is given m 
operator as 


34| . The h-function interaction can be expressed through the density 


Sint Cl 


N1N2 

N 


dttr(0i02), 


(58) 


where a = g/y/N, and 0* {i = 1,2) correspond to different spin directions. Here we are 
taking the limit that for each spin direction there are a large number of particles, Ni and N 2 
respectively, and both of them are of the same order as the total number Y {q]. In writing 
the above formula we have implicitly used the classical nature of the density operator in the 
large-Y limit, 

{x\^\y) = (j){x)6{x-y), (59) 


which is simply the coordinate representation of fl3Tl) . One can check that if (pi and 02 are 
identical, the interaction is indeed trivial due to dsa and (1501) . Now with the interaction, 
the equations for 0i and 02 become entangled 

idt^i + [h, 0i] = a^^^[02, 0i] 

idt^2 + [h, 02] = a ^ [01, 02]- (60) 

In the hrst equation the back-reaction of 02 to 0i is proportional to the hnite ratio Y 2 /Y. 
The Yi-factor cancels out since the corresponding Fermi energy scales as Yi. The conclusion 
is also true for the second equation. Therefore only the hnite ratios remain in the large-Y 
limit, as one expects. 

In the limit a —?• 0 they decouples, and one recovers two free-fermion systems. For small a 
the coupled equations could be solved perturbatively. For an inhnitely repulsive interaction, 
minimizing the action (1581) forces the product of the two density operators to vanish 


0102 — 0 . 


(61) 
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In other words, fermions of different spin direction can not occupy the same energy 
state/phase space region any longer. This gives rise to an effective exclusive condition be¬ 


tween different spin components 


21 


33[ |. Therefore the total density operator, 0 = < 


describes simply a free-fermion system with no distinguishing of the spin occupation. That 
means, 0 satishes the constraints 05111501) and the free equation fl5TD . and the solution in 
phase space is given by flFH) . Notice that such a conclusion is actually independent of the 
exact form the single-particle Hamiltonian ([53]). The situation for general g will be com¬ 
plicated. Using the phase space density (l){p,q), the equations fl60D can be reexpressed in 
terms of the Moyal bracket, which reduces to the Poisson bracket in the large-limit. Since 
the Poisson bracket involves only double derivative, the equations seem to be tractable. It 
would be interesting to see if in this way one can recover some kind of generalized Fredholm 
equation 2, 


V. SUMMARY 

In the paper we have shown how to obtain the planar limit of one dimensional many-body 
system. For the free fermion system it inherits from the matrix theory. In the bosonic case, 
the interactions can be introduced as the potential term of the density held, the collective 
held. Such a hied behaves classically in the planar limit, and can be analyzed easily. For the 
fermions, only interactions between diherent components survive. In the planar limit they 
can be expressed as the overlap integral of the phase density functions. Evolution of the 
system is determined by the corresponding classical equations in terms of Poisson brackets. 

We use the ^-function interaction to show how these methods can be explicitly applied. 
Such a model exhibits a dynamical evolution from an interacting system to a free fermionic 
one. This provides a dynamical bosonization of nonrelativistic fermions, and shows the 
diherences between the collection held theory and the fermion phase description. Without 
using the Fredholm equation from the Bethe ansatz, we recover the correct results at strong 
and weak interacting limits. In this sense, the present formalism serves as a general com¬ 
plementary description to the Bethe ansatz, which is specihc for the model. In particular, 
the planar limit gives a solid basis for the validation of the localized Fredholm equation 
proposed recently. The methods could be immediately applied to general interactions, and 
even higher dimensional systems. Also one could include the sub-leading terms to study the 
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finite-corrections. One may expect such corrections give rise to the low energy fluctua¬ 
tions around the ground state density. We will try to investigate some of these topics in the 
future. 
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